function  N=f_u(u);
N=3.*u.^2;